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We present the most accurate and complete data set for the analyzing power A y (6) in neutron- 
proton scattering. The experimental data were corrected for the effects of multiple scattering, both 
in the center detector and in the neutron detectors. The final data at E n = 12.0 MeV deviate 
considerably from the predictions of nucleon-nucleon phase-shift analyses and potential models. 
The impact of the new data on the value of the charged pion-nucleon coupling constant is discussed 
in a model study. 

PACS numbers: 13.75.Cs,24.70.+s,25.40.Dn 



INTRODUCTION 

There are reasons why low-energy (E^ < 20 MeV) 
Nucleon-Nucleon (NN) scattering data might appear to 
be of limited use in constraining NN phase-shift analyses 
(PSAs) and potential models (PMs) 0, 11 @, S> • 

For one thing, the deuteron bound-state properties al- 
ready provide a fairly stringent constraint for any NN 
PM, and might seem sufficient. For another, low-energy 
scattering data can provide constraints only for the lower 
partial wave NN interactions and, even then, can not 
determine individual partial waves. For example, the 
low-energy analyzing power, A y (8), is governed by the 
three angular momentum L = 1 interactions, 3 Po, 3 -Pi, 
and 3 P2- Although NN data provide constraints on the 
3 P phase shifts taken together, it cannot determine each 
parameter unambiguously Q. 

Despite the very small magnitude of NN A y (0), its 
importance derives from the fact that it is possible ex- 
perimentally to measure such data to great precision. As 
a result, NN A y {&) can provide a crucial test of our un- 
derstanding of the NN interaction and of the nuclear 
force in general. 

Perhaps the most important controversy surround- 
ing current NN interaction models concerns the pion- 
nucleon coupling constant, g^/Air. In Ref. [9], Machlcidt 
and Slaus point out that low-energy proton-proton (p-p) 
A y (8) data are very sensitive to the neutral pion-nucleon 
coupling constant, implying a value of g^o/47r < 13.4 (see 
also Ref. [10|). At the same time, the correct descrip- 
tion of the quadrupole moment of the deuteron and low- 
energy neutron-proton (n-p) A y (9) data requires meson- 
exchange based TV TV potential models to have values for 
the neutral and charged pion-nucleon coupling constants 
<7^o/47r and s£±/47r, respectively, of 14.0 or larger. The 
latter finding is clearly inconsistent with the results of 
the Nijmegen group's NI93 PSA (^ /4tt = 13.47 ± 0.11 



and gl ± /4n = 13.54 ± 0.05 11]) and of the VPI group 
{rLo/±k = 13.3 and s£ ± /4tt = 13.9 from NN scattering 
JJ, til] and 5 2 /4tt = 13.75 ± 0.15 from w ± p scattering 

m 



In principle, this inconsistency can be reduced by as- 
suming a charge-splitting of the pion-nucleon coupling 
constant, i.e., one could assume that the neutral pion 
couples to the nucleon with a different strength than the 
charged pion. In Refs. [9j, [l0| it was shown that the com- 
bination of g 2 /47T = 13.6 and g^±/47r = 14.4 creates a 
sufficiently large value for the quadrupole moment of the 
deuteron, and reproduces the low-energy p-p 3 Pq phase 
shifts. 

At the same time that the analyses of Refs. 0, [l(| 
were performed, there were indications that n-p differen- 
tial cross-section data at intermediate energies favored 
a larger value for <7^ ± /47r. On this, see Ref. [IH for 
a comprehensive overview of recent determinations of 
5^/47r and especially Ref. [3], which quotes g^/4-ir = 
14.50 ± 0.26 obtained from n-p differential cross section 
data at E n = 162 MeV. However, the recent n-p differ- 
ential cross section data obtained at IUCF at 194 MeV 
do not support this larger value of g^ ± /47r. 

Although it seems likely that there is no significant 
charge splitting in the pion-nucleon coupling constants 
at intermediate energies, the question remains unresolved 
at low energies. On the one hand, the theoretical models 
used to account for the charge dependence of the singlet 
NN scattering lengths, 1 Sq, do not allow for any large 
charge splitting of g^/4ir. On the other hand, many low- 
energy data suggest a significant charge splitting. We 
report here on the results of a new n-p A y (8) experiment 
carried out at E n — 12.0 MeV utilizing improved data- 
taking and data- analysis techniques. For references to 
previous n-p A y {8) measurements see Refs. 3, Ijl 20| . 
Our results confirm the inconsistencies between low- 
energy analyzing power and available theoretical models 
of the NN interaction. 



EXPERIMENTAL SETUP 

The experimental setup is shown in Fig. [1] Polarized 
neutrons with mean energy of 12.0 MeV and total en- 
ergy spread of about 400 keV were produced via the 
polarization-transfer reaction 2 H(d,n) 3 He at 0°. The po- 
larized deuteron beam was accelerated to Ed = 9.40 MeV 
and entered through a 4.6 /im Havar foil into a small (3.14 
cm long, 0.48 cm radius) gas cell filled with 7.8 atm of 
deuterium gas and capped with a 0.1 cm thick gold beam- 
stop. The gas cell was mounted inside a 1.8 m-thick wall 
made of concrete, paraffin, iron, copper, and lead, to 
shield the neutron detectors from the direct flux of the 
neutron source. Typical deuteron beam intensities on 
target were 1.5 /iA and typical values for the deuteron 
vector polarization \p z \ were 0.65. The polarized neu- 
trons produced at 0° relative to the incident deuteron 
beam passed through a collimation system to produce a 
rectangular shaped neutron beam at the position of the 
proton-containing active target labelled "Center Detec- 
tor" in Fig. [TJ The center detector (CD) consisted of 
an upright cylinder made of the plastic scintillator ma- 
terial NE102A with dimensions 1.9 cm diameter and 3.8 
cm height. The CD was located at a distance of 172 cm 
from the neutron source and was mounted via a short 
light guide onto a 5 cm diameter photomultiplicr tube 
(PMT). 

Neutrons scattered to the left or right were detected 
by five pairs of neutron detectors positioned symmetri- 
cally relative to the incident neutron beam direction in 
the horizontal scattering plane. The neutron detectors 
(NDs) were filled with the liquid scintillator material 
NE213. These detectors had excellent neutron-gamma 
pulse-shape discrimination capabilities and had an ac- 
tive volume of 4.3 cm wide, 11.9 cm high and 7.5 cm 
deep. They were viewed by 5 cm diameter PMTs through 
0.5 cm thick Pyrex glass windows and 7.5 cm long light 
guides. The neutron detectors were mounted onto (low- 
mass) 30 cm high stands and placed on an aluminum ring 
surrounding the CD. The center-to-center distance be- 
tween the CD and the neutron detectors ranged from 45 
cm to 70 cm depending on scattering angle. The angular 
separation between the neutron detectors was 12° (lab). 
In order to cover the angular range from 9i a b = 16° to 72° 
in 4° steps, three settings of the five detector pairs were 
required. The absolute magnitude of the neutron polar- 
ization was measured with a neutron polarimeter located 
downstream of the n-p scattering arrangement. The po- 
larimeter consisted of a 4 He gas scintillator pressurized to 
100 atm (95% He, 5% Xe) and a pair of neutron detectors 
positioned at 6i a b — 58°, which were identical to those 
used for n-p scattering. In order to reduce instrumental 
asymmetries for the n-p and n- 4 He measurements, the 
deuteron vector polarization p z , and therefore the neu- 
tron polarization, was flipped at a frequency of 10 Hz 



(between up and down relative to the horizontal scatter- 
ing plane). The n-p and n- 4 He data were accumulated 
simultaneously in six runs, each lasting about 250 data- 
taking hours. 

The data-acquisition electronics recorded the center 
detector pulse height (CDPH), the neutron time of flight 
(NTOF) between the CD and the NDs, and spectra for 
each ND, including pulse-shape information. Since the 
energy of the scattered neutrons varied from E n i = 11.1 
MeV at 6 lab = 16° to E n , = 1.1 MeV at B lah = 72°, 
different hardware thresholds were used for the NDs. In 
addition, three different gains were used for the CD sig- 
nals (using different dynodes). 

Software cuts were set on the CDPH and the pulse 
height in the NDs to eliminate pulses at the extreme ends 
of the spectra. Gates were also set on the neutrons in the 
pulse-shape discrimination spectra and wide gates were 
set on the elastic peak of the NTOF spectrum. All four 
of these cuts (identical for spin-up and spin-down spec- 
tra) were used to generate two-dimensional (2D) spectra 
of CDPH versus NTOF, for scattering to the left and 
right NDs and for neutron spin up and spin down. An 
example of such a spectrum is shown in Fig. [2] where the 
CDPH scale has been temporarily compressed in order 
to fit within 64 channels. Tight NTOF gates were set 
in these 2D spectra eliminating the tails of the peak, as 
shown in Fig. [2 in order to identify the elastic scattering 
events of interest (again, identical for spin-up and spin- 
down spectra) . These new NTOF gates were used to sort 
the final CDPH spectra (now in their full 512-channel res- 
olution) corresponding to each neutron detector and spin 
state. The CDPH spectra were used to determine the n- 
p yields and scattering asymmetries, after applying the 
corrections described in the following section. The above 
process was also followed to sample the accidental (i.e., 
time uncorrelated) background by using an NTOF cut 
located at times shorter than the gamma peak. The ac- 
cidental background proved to be extremely small. 

DATA ANALYSIS 

After the sorting procedure described above and the 
subtraction of the accidental events, the data still 
contained a number of finite-geometry and multiple- 
scattering effects. To remove these effects, Monte-Carlo 
calculations were performed to simulate the experiment. 
Two effects are due exclusively to the finite size of the 
center detector (CD) and the neutron detectors (NDs) 
and have a slight effect on single-scattering events. First, 
because there is a range of angles subtended by each de- 
tector set at each nominal angle and because the cross 
section of n-p elastic scattering varies over this range, we 
must report an effective angle. These were calculated by 
our code and are listed in the first column of Table [J 
This effect is small; the largest shift is no more than a 



half of a degree. The second finite-geometry effect con- 
cerns the value of A y (9) itself, again due to the range 
of angles subtended by each ND. Effective A y {9) values 
were calculated by our Monte-Carlo code and these were 
compared to the values from the code's library. The ratio 
between these two values was then applied to the data. 
Once again, the correction is small; only the first four an- 
gles had corrections that were larger than the uncertainty 
of the calculation (about 0.00012). 

In addition to elastic scattering, multiple scattering 
events occur in the CD. About 50% (depending on ND 
angle) of these events were eliminated as a result of the 
neutron time of flight (NTOF) gate. Nevertheless, the 
center detector pulse height (CDPH) spectrum contained 
multiple scattering events amounting to approximately 
2% of all single scattering events. Our Monte-Carlo sim- 
ulation showed that the only significant processes were 
those due to double scattering, specifically neutron dou- 
ble scattering from hydrogen ( 1 H- 1 H), neutron scatter- 
ing from hydrogen and subsequent scattering from car- 
bon ( 1 H- 12 C), and neutron scattering from carbon and 
subsequent scattering from hydrogen ( 12 C- 1 H). In per- 
forming these calculations, we used complete libraries of 
cross-section and polarization data for both n- 1 !! and n- 
12 C scattering. We will return to the subject of the n- 12 C 
library in our discussion of the PDE correction. 

We also removed edge-effect events from the data, 
which result when recoil protons leave the CD before de- 
positing their full energy. Along with the double scatter- 
ing events, these counts elongate the tails of the CDPH 
peak, especially to the left (low-energy) side. 

A sample CDPH spectrum is shown in Fig. [3] (top 
panel). The solid curve represents our Monte-Carlo sim- 
ulation for a scattering angle of #; a & = 36°, while the 
small open circles show the experimental data. A greatly 
expanded view is shown in the middle panel, where the 
open circles again indicate the experimental data. The 
curves labeled "double" are the calculated double scatter- 
ing contributions 1 H- 1 H, l H- 12 C, and ^C^H. The dotted 
curve labled "edge" is the calculated pulse-height distri- 
bution due to edge effects. Finally, the curve labeled 
"single" is the calculated single scattering contribution 
(plus the edge effects). All of the calculations are nor- 
malized to the data. 

Another expanded view is displayed in the bottom 
panel of Fig. [3l Here, "data with subtraction" shows the 
data after the removal of all counts due to double scat- 
tering and edge effects (labeled "ms+edge" ) . Even after 
this subtraction, a small background remains, amounting 
to about 0.3% of the single-scattering events. A number 
of fits were used to estimate this remaining background, 
ranging from a linear fit between channel numbers 150 
and 350 to a parabolic fit between channel numbers 180 
and 280. Due to the smallness of the remaining back- 
ground, the asymmetry proved to be independent of our 
background choice, within statistical uncertainties. We 



also concluded that the background was unpolarized. For 
all ND angle settings we approximated the remaining 
background by a linear function connecting the left and 
right sides of the CDPH peak (for example, in Fig.[3]from 
channel 180 to 280). The remaining background seen 
above channel 290 in the bottom panel of Fig.[3]is due to 
cross-talk effects between two adjacent detectors, specif- 
ically neutron scattering from the detector positioned at 
a larger scattering angle (and shorter distance from the 
CD) to the detector of interest. 

Three sets of gates were used to calculate the yields 
and asymmetries, at 10% (shown in Fig. [3]by the dashed 
lines), 30%, and 50% of the CDPH peak maximum. In 
Fig. [31 this is done for , for spin up scattering to the 
left ND at 9 iab = 36°. Similarly, the yields N^, N^, 
and Nft were obtained to calculate the asymmetry e = 

(a — l)/(a + 1) with a = \/ttt"^t- The nominal gates 

y "r "l 

were the 30% set. The other two gates (the 10% and 
50% set) were used to check on the appropriateness of the 
background subtraction. Within statistical uncertainty, 
the results for e proved to be independent of the choice 
of the gate width. 

In order to extract the n-p A y (9) from the measured 
asymmetry e(9), the neutron polarization p y must be 
known. For this purpose the n- 4 He asymmetry data ac- 
quired with the neutron polarimeter referred to above 
were processed and analyzed in the same way as the 
n-p asymmetry data. In this case the 4 He recoil pulse 
height in the high-pressure gas scintillator plays the role 
of the CDPH in the plastic scintillator used for the n- 
p asymmetry measurements. The neutron polarization 
was obtained from eij e (58°) = {otHe — l)/(<XHe + 1) — 
^(58°)^™, where an e is defined as above. Here, the ef- 
fective analyzing power A y (58°) for n- 4 He scattering at 
E n = 12.0 MeV was calculated for the present neutron 
polarimeter geometry via Monte-Carlo calculations. The 
n- 4 He phase shifts of Stammbach and Walter [2l| were 
used. All of the relevant multiple scattering processes 
were included. We obtained A y {58°) = -0.554 ± 0.008, 
where the uncertainty is mainly of a systematic na- 
ture reflecting the uncertainty associated with the n- 
4 He phase shifts. The average neutron polarization was 
p n y = 0.563 ± 0.008. 

At such a high level of precision, a subtle systematic ef- 
fect comes into play, which does not cancel by reversal of 
the neutron polarization. This is the polarization depen- 
dent efficiency (PDE) [n| of the neutron detectors. The 
NDs contain hydrogen and carbon in the ratio of 1.21:1. 
The double scattering process 12 C- 1 H in the NDs, which 
accounts for about 10% of the total neutron detection 
efficiency, is sensitive to the n- 12 C A y (6). If the n- 12 C 
A y (9) is not constant over the range of neutron energies 
E n i seen by a particular ND, an instrumental asymmetry 
will occur. Typical values for AE n r are 800 keV. A realis- 



tic correction for this effect requires a detailed knowledge 
of the n- 12 C A y (8), especially in the resonance region of 
the n- 12 C total cross section between 2.0 and 8.5 MeV 
neutron energy. In this energy regime the n- 12 C A y (9) 
changes rapidly and therefore causes sizeable PDE ef- 
fects. 

All of the post-1985 n-p A y (0) measurements have been 
corrected for the PDE. However, due to the lack of a de- 
tailed n- 12 C A y (8) database, especially at low energies, 
the accuracy of the associated corrections was limited. 
In assembling our data library, we used the thirty-three 
n- 12 C A y (9) angular distributions measured by Roper et 
al. [13] in the energy range from 2.2 to 8.5 MeV. From 
E n = to 6.5 MeV, we used an R-matrix analysis by Hale 
[23( 1 . which included the data from Ref. [22j]. In certain 
regions (especially for forward angles and for neutron en- 
ergies between 3.5 and 4.5 MeV), the analysis of Ref. (23| 
missed the A y {9) data slightly and we therefore substi- 
tuted Legendre polynomial fits to the data of Ref. [22| in 
these regions. Between 6.5 MeV and 8.5 MeV, we used 
fits to the data of Ref. [22| as well as the recent phase- 
shift analysis (PSA) of Chen and Tornow 24(. Above 
E n = 8.5 MeV, we used the Chen- Tornow PSA exclu- 
sively. The new data by Roper et al. and the analyses 
of Hale, Chen and Tornow improved the n- 12 C A y {9) 
database considerably, making corrections for the PDE 
more reliable. 

We ran our Monte-Carlo code for 20 separate legs, each 
leg of three million events, and each leg starting from a 
different random number. The PDE correction to the 
A y (6) data was taken as the difference between the A y (9) 
result with polarization effects turned on in the NDs and 
the result with the polarization turned off. The second 
column of Table |l] lists our final PDE corrections. Note 
that they vary greatly from one data point to the other 
due to the pronounced resonance features in n- 12 C scat- 
tering at low energy. It is important also to note that 
our present results agree well with the overall trend of 
the PDE corrections of Ref. 20] , which used a different 
Monte-Carlo code and a different database. Our reason 
for having much greater confidence in the present PDE 
results is due to our extensive and detailed work in re- 
vising the data libraries, as outlined above. 

The third column in Table H] summarizes our final re- 
sults for A y (6) in n-p scattering at E n = 12.0 MeV. Note 
the small overall uncertainty. The final results include 
uncertainties in A y {ff) due to statistics, the measure- 
ment of beam polarization, the multiple-scattering cal- 
culation, the PDE calculation, and the remaining back- 
ground (typically zero) all added in quadrature. The 
final uncertainties are about half of those of the previous 
TUNL n-p A y {0) measurement at E n = 12.0 MeV 
This is partly due to the fact that the Atomic Beam Po- 
larized Ion Source used in the present study produced 
about four times the deuteron current as the Lamb-Shift 
source used in the previous study. 



DISCUSSION 

Figure f?] shows the present n-p A y (9) in comparison 
to the NN phase-shift analysis prediction (solid curve) 
of the Nijmegen group, NI93. Clearly, NI93 provides a 
larger A y {9) throughout the entire angular distribution. 
The accuracy of the neutron polarization determined in 
the present work does not allow for a renormalization 
of the A y {9) data beyond the error bars given in Fig. 
[4] Furthermore, the present n-p A y (9) data are in good 
agreement with the trend established by previous TUNL 
data where a different method was used for determining 

As we have pointed out in the introduction, the un- 
derlying NN dynamics that characterize A y (9) precludes 
us from extracting unambiguous information about the 
3 Pj NN interactions. However, we can conclude that the 
NI93 NN PSA overestimates the n-p A y (9) at E„ = 12.0 
MeV. This statement is of considerable importance con- 
sidering the fact that most NN potential model builders 
use the NI93 PSA results or the associated database for 
determining the free parameters of their models. One has 
to conclude that all the recent so-called high-precision 
NN potential models overestimate the n-p A y (9) at low 
energies. This observation has far-reaching consequences 
for nuclear scattering systems with A > 2, which are 
much more sensitive to the 3 Pj NN interactions than 
the NN system [Hj]. 

Valuable information can be obtained from the present 
data if they are compared to variations of the theoretical 
predictions. Here we focus on the charged pion coupling 
constant [^, [l(| ■ Figure 0] shows our data in compar- 
ison to three theoretical predictions based on the CD- 
Bonn NN potential, which use three different values of 
the charged pion-nucleon coupling constant, g 2 ± /4ir. In 
these three models, only the S-wave NN interactions of 
CD-Bonn were refitted. All three predictions use the 
same neutral pion coupling constant, g 2 /47r = 13.6. 
The curve using g 2 ± /4ir = 13.6 is indistinguishable on 
this scale from the prediction of NI93 (solid curve) . The 
dashed curve in Fig.|4]uses g 2 ± /Air = 14.0 and the dotted 
curve uses g 2 ± /An = 14.4. The values of x 2 per degree 
of freedom associated with the solid, dashed and dotted 
curves are 6.0, 1.7, and 2.5, respectively. Therefore, this 
model study confirms and puts on more solid ground the 
findings of Refs. [^, [1(3] regarding low-energy n-p A y {9) 
data and their demand for a charge splitting of g 2 /An. 

In summary, the present data represent the most ac- 
curate and complete n-p A y {9) angular distribution ever 
reported. Our model study based on the CD-Bonn NN 
potential model supports a substantial charge depen- 
dence of the pion-nucleon coupling constant. Our re- 
sults are inconsistent with the existing global NN PSAs 
of the Nijmegen [ll[ and VPI II, 13] groups and with 
high-precision NN potential models. However, our re- 



suits agree with inconsistencies previously noticed be- 
tween data and predictions for the 3 S\- 3 Di mixing pa- 
rameter t\ in n-p scattering at low energies [261 ] and also 
with requirements placed on the charged coupling con- 
stant by the quadrupole moment of the deuteron [9( . Of 
course, it is possible that neither of these scenarios is 
the "correct" one. Perhaps the impasse comes because 
we are at the point where the precision of our data and 
the development of our "low energy" theoretical models 
has pushed the paradigm of meson-exchange based NN 
potential models beyond its limits. 

This work was supported in part by the U.S. Depart- 
ment of Energy, Office of Nuclear Physics, under Grant 
No. DE-FG02-97ER41033. The authors would like to 
thank R. Machlcidt for valuable contributions to this 
work. 



TABLE I: Results of n-p A y (8) experiment at E n = 12.0 MeV. 





PDE correction 


final results 


32.6 


0.00014 ± 0.00016 


0.00854 


± 


0.00067 


40.5 


0.00004 ± 0.00016 


0.01231 


± 


0.00064 


48.5 


0.00006 ± 0.00015 


0.01451 


± 


0.00065 


56.5 


0.00013 ± 0.00013 


0.01443 


± 


0.00063 


64.4 


-0.00005 ± 0.00015 


0.01560 


± 


0.00063 


72.4 


0.00294 ± 0.00022 


0.01659 


± 


0.00067 


80.5 


-0.00185 ± 0.00014 


0.01470 


± 


0.00060 


88.4 


0.00019 ± 0.00017 


0.01386 


± 


0.00057 


96.3 


0.00072 ± 0.00018 


0.01198 


± 


0.00059 


104.2 


-0.00136 ± 0.00018 


0.01110 


± 


0.00058 


112.2 


-0.00114 ± 0.00028 


0.00662 


± 


0.00062 


120.2 


0.00108 ± 0.00029 


0.00558 


± 


0.00065 


128.2 


0.00103 ± 0.00021 


0.00483 


± 


0.00056 


136.0 


-0.00018 ± 0.00036 


0.00372 


± 


0.00067 


143.8 


-0.00029 ± 0.00040 


0.00287 


± 


0.00079 
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FIG. 1: Experimental setup for n-p A y (8) measurements in 
TUNL's Shielded Neutron Source area. 





FIG. 2: 2D spectrum of compressed CDPH versus NTOF 
for scattering to #; a j, = 64°. A tight NTOF gate was set 
around the elastic neutron peak in order to remove as many 
background events as possible. 
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FIG. 4: Neutron-proton A y {d) data at £ n = 12.0 MeV in 
comparison to theoretical predictions. The error bars asso- 
ciated with the data represent the overall uncertainty of the 
data with statistical and systematic uncertainties added in 
quadrature. The solid curve is the Nijmegen NI93 PSA pre- 
diction. The other curves are for the CD-Bonn based model 
study which varies the charged pion coupling constant. Here, 
for <?^o/47r, all three curves use 13.6. For </^±/47r, the cal- 
culation using 13.6 coincides on this scale with the Nijmegen 
NI93 PSA result (solid curve); the dashed curve uses 14.0 and 
the dotted curve 14.4. 



FIG. 3: The top panel shows a comparison of calculated (solid 
curve) and measured (dots) CDPH spectrum for scattering to 
Oiat = 36°. The middle panel shows an expanded view with 
focus on calculated multiple-scattering and edge effect con- 
tributions. The bottom panel shows an additional expanded 
view focusing on the remaining background. See text for de- 
tails. 



